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^  CONGRATS  (CONtinuous  Gradient  RAy  Tracing  System)  is  an  integrated 
collection  of  ray  tracing  programs  designed  to  model  acoustic  wopaga- 
tion  and  reverberation^'lis  described  in  references  (&)  and  (b)7^Al- 
though  sound  speed  data  is  often  obtained  from  temperature  and  salinity 
readings  at  various  depths,  the  fundamental  CONGRATS  programs,  S0990 
and  SO9910  formerly  required  a  table  of  velocity  versus  depth  as  an 
input.  Hence,  it  was  sometimes  necessary  to  convert  the  empirical  data 
into  a  velocity-depth  profile  before  using  the  CONGRATS  series.  CONGRATS 
has  now  been  revised  to  convert  a  temperature  and  salinity  profile  to_a 
velocity  profile  automatically^  Velocity  data^CTn  the  format  described 
previously  in  reference  (a),''^2B  still  accepted  by  the  programs.  It  is 
hoped,  however,  that  the  use  of  the  new  option  will  remove  the  burden 
of  some  preliminary  hand  computations  from  the  user.  This  memorandum 
-w±3±  containAh  discussion  of  the  metnod  used  in  the  conversion  of  tem¬ 
perature  to  velocity,  a  detailed  description  of  the  data  necessary  to 
implement  the  conversion,  a  sacgjle  run  along  with  its  output,  and  a 
listing  of  the  revised  Subroutine  INPUT  and  the  new  Subroutine  BT.  a 


METHOD  ' 

Although  Wilson's  equation  is  most  widely  used  in  the  calculation 
of  sound  speed  in  vrater,  a  sing)lified  formula  was  programmed  for  the 
CONGRATS  series.  According  to  C.  C.  Leroy,  this  formula  fits  Wilson's 
data  with  a  better  accuracy  than  does  V/ileon's  equation  (over  a  domain 
restricted  to  areas  of  operational  interest)  and  approaches  Greenspan 
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and  Tschiegg's  results  for  very  low  salinity  water  (see  reference  (c)). 

The  fonmila  is  presented  in  Table  1. 

The  basic  formula,  V  =  ,  is  sufficiently  accurate  (see 

reference  (c))  tinder  the  following  conditions:  depth  less  than  "J ,000 
meters,  temperature  less  than  25°  Centigrade,  and  salinity  between  30 
and  40  parts  per  thousand.  The  remaining  three  terms  are  added  for 
increased  accuracy  if  the  above  conditions  are  not  met.  If  the  tem¬ 
perature  is  greater  than  or  equal  to  25°  centigrade,  the  correction 
term  is  added;  if  the  depth  is  greater  than  7 >000  meters,  is 

added;  and  if  the  salinity  is  less  than  30  parts  per  thousand,  Vd  is 
added. 

DATA  DESCRIPTION 

The  input  deck  of  CONGRATS  Programs  SO99O  and  SO99I,  as  described 
in  reference  (a),  consists  of  sets,  each  of  which  is  an  ordered  collec¬ 
tion  of  standard,  eighty- column,  punched  cards.  The  conversion  of  tem¬ 
perature  and  salinity  data  to  a  velocity  profile  is  implemented  by 
substituting  a  ten^ierature  and  salinity  profile  set  for  the  velocity 
profile  set.  The  rest  of  the  input  deck  is  unaffected  by  this  change. 

The  first  card  of  the  temperature  and  salinity  profile  set  is  divided 
into  six  fields  of  ten  columns  each,  following  the  format  convention 
set  in  reference  (a).  The  first  field  contains  the  word  "THERMAL," 
starting  in  column  1,  and  the  second  contains  the  word  "PROFILE," 
starting  in  colvimn  11.  These  two  fields  identify  the  temperature  and 
salinity  profile  set.  The  third  field  contains  either  the  word 
"CONTINUOUS,"  in  which  case  the  resultant  velocity  profile  is  fitted 
with  continuous  gradients,  or  the  word  "CONSTANT,"  in  which  case  the 
constant  gradient  curve  fitting  technique  is  used.  In  either  case, 
the  word  mu.st  begin  in  column  21.  The  last  three  fields  of  the 
THERI»1AL  PROFILE  card  are  numeric  fields,  each  of  which  is  read  into 
the  computer  using  an  FIO.5  format.  The  first  numeric  field  contains 
the  number  of  points  in  the  profile  and  the  second  numeric  field  con¬ 
tains  the  latitude  in  degrees.  The  third  numeric  field  may  contain 
the  salinity  in  parts  per  thousand,  if  the  salinity  is  constant  through¬ 
out  the  profile.  If  salinity  varies  with  depth,  this  field  may  be  left 
blank. 

The  second  card  of  the  ten5)erature  and  salinity  profile  set  is  the 
units  card  which  uses  three  fields  of  ten  columns  each,  starting  in 
columns  1,  11  and  21  respectively.  The  first  field  contains  the  units 
of  depth,  the  second  contains  the  units  of  temperature  and  the  third 
indicates  the  units  in  which  the  resultant  velocities  are  to  be  printed  ^4':^ 
and/or  plotted.  Table  2  of  reference  (a),  together  with  the  following 
additions,  is  a  complete  list  of  acceptable  input  units,  their  CONGRATS 
abbreviations,  their  conversion  factors  and  the  resulting  program  units. 
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Input  Units 

Abbreviation 

Conversion  Factor 

Result 

Centigrade 

C 

i.ooooopoo 

C 

Fahrenheit 

Fahr 

0.55555555 

C 

Note:  32  must  be 

subtracted  from  the 

number  of  degrees  Fahrenheit 

before  multiplication  by  the  conversion  factor. 


The  cards  containing  the  temperature  profile,  arranged  in  ordar 
of  increasing  depth,  follow  the  units  card.  These  cards  use  three 
numeric  fields,  starting  in  columns  1,  11  and  21,  of  ten  columns  each: 
the  first  field  contains  the  depth,  the  second  contains  the  teii5)erature 
and  the  third  contains  the  salinity.  If  the  salinity  is  constant  as 
a  function  of  depth  and  the  salinity  field  has  been  filled  on  the  first 
card  of  the  set,  the  third  field  of  the  profile  cards  may  be  left  blank 

The  possible  choices  for  Indicating  the  salinity  are  charted  in 
Table  2.  If  the  seilinity  field  on  each  card  of  the  temperature  profile 
and  the  constant  salinity  field  on  the  first  card  of  the  set  are  both 
left  blank,  then  the  salinity  is  considered  to  be  zero  parts  per 
thousand  throvighout  the  profile.  If  the  salinity  field  on  each  profile 
card  is  left  blank  but  the  constant  salinity  field  on  the  first  card  of 
the  set  contains  a  positive  nvunber,  then  the  salinity  is  considered 
constant  and  its  value  is  the  number  indicated  on  the  first  card.  In 
both  cases  the  resultant  constant  salinity  value  is  printed  at  each 
depth  of  the  temperature  and  salinity  profile  in  the  computer  print¬ 
out.  If  a  positive  value  is  shown  in  the  third  field  of  one  or  more 
profile  cards,  the  constant  salinity  field  of  the  first  card  is  ignored 
and  the  resultant  salinity  profile  consists  entirely  of  the  salinity 
values  from  the  profile  cards.  It  should  be  noted  that  a  negative 
salinity  will  cause  the  program  to  terminate  with  an  error  stop. 

The  cards  presented  in  Fig.  1  are  an  example  of  a  temperature  and 
salinity  profile  set.  The  first  card  indicates  that  the  resultant 
velocity  profile  is  to  be  fitted  with  continuous  gradients,  that  there 
are  I9  points  in  the  profile,  that  the  latitude  is  39  degrees,  and  that 
the  salinity,  which  is  constant,  is  38*2  parts  per  thousand.  The 
second  card  indicates  that  depths  are  in  feet,  ten^jeratures  are  in 
degrees  Fahrenheit,  and  that  the  resultant  velocities  are  to  be  printed 
in  feet  per  second.  The  remaining  cards  contain  the  depths  and  tem¬ 
peratures  of  the  profile.  Sedinity  values  are  absent  from  these  cards 
because  the  salinity  is  constant  and  indicated  on  the  first  card  of  the 
set. 


An  example  of  a  temperatxire  and  salinity  profile  in  which  salinity 
varies  with  depth  is  presented  in  Fig.  2.  The  use  of  the  continuous 
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gradient  curve  fitting  technique  is  requested;  there  are  *10  points  in 
the  profile,  and  the  latitude  is  40°.  Bae  depth  units  are  feet,  the 
temperatui’e  units  are  degrees  Fahrenheit,  and  the  velocity  units  are 
feet  per  second.  Since  salinity  varies  with  depth,  the  profile  cards 
indicate  temperature  and  salinity  for  each  depth  and  the  salinity  field 
on  the  first  card  of  the  set  is  left  blank. 

If  a  plot  of  temperatvire  versus  depth  and/or  salinity  versus  depth 
is  desired,  a  THERMAL  AXES  card  must  be  included  in  the  input  deck. 

Hiis  cenTd  contains  the  words  "THERMAL"  and  "AXES"  starting  in  columns 
1  and  11,  respectively.  The  third  field  contains  the  units  in  which 
the  axes  are  to  be  plotted.  Inches  and  centimeters  are  the  available 
units.  The  fourth  field,  which  is  numeric,  contains  the  length  of  the 
depth  axis,  which  is  plotted  vertically.  The  numbers  in  the  fifth  and 
sixth  fields  are  the  lengths  of  the  horizontal  temijerature  and  salinity 
axes,  respectively.  A  zero  ten^jerature  axis  length  or  salinity  axis 
length  will  sujjpress  the  temperat\a^  versus  depth  or  salinity  versus 
depth  plot,  respectively.  Both  plots  will  be  suppressed  if  the  depth 
axis  length  is  zero.  Figure  3(a)  presents  a  THERMAL  AXES  card  which 
would  cause  both  a  temx)erature  and  salini\;y  profile  to  be  plotted.  The 
depth  axis  would  be  10  Inches  long  in  both  plots;  the  temperature  axis, 

8  inches  long;  and  the  salinity  axis,  5  inches  long.  The  card  shown 
in  Fig.  3(b)  wovild  cause  only  one  plot  to  be  drawn.  The  temperature 
versus  depth  plot  would  have  a  depth  axis  of  10  inches  and  a  temperature 
axis  of  6  inches,  and  the  value  of  the  salinity  at  the  surface  would  be 
printed  at  the  left  of  the  plot.  If  the  user  desires  other  information 
plotted  by  the  program,  the  appropriate  AXES  cards  (as  described  in 
reference  (a))  must  be  added  for  each  type  of  plot,  e.g.,  for  a  velocity 
profile  plot  to  be  drawn,  a  VELOCITY  AXES  card  must  be  added  to  the 
input  deck. 


EXAMPLES 

Two  different  sample  runs  have  been  selected  to  illustrate  the  use 
of  a  teraperatm'e  and  salinity  profile  in  the  CONGRATS  series.  Example  1 
uses  the  temperature  and  salinity  profile  shown  in  Fig.  1.  Figure  4  is 
a  listing  of  the  Example  1  run  deck.  The  resultant  computer  print¬ 
out  (see  Fig.  5)  shows  the  temperature  and  salinity  profile  with  the 
constant  salinity  value,  38*2  parts  per  thousand,  printed  along  with 
each  depth  and  temperature  of  the  profile.  The  latitude,  which  is  39°, 
is  printed  below  the  thermal  profile.  The  velocity  profile,  which  has 
been  computed  by  the  program,  is  listed,  followed  by  the  velocity 
tolerance  used  to  fit  the  data.  Tie  THERMAL  AXES  card  (see  Fig.  3(b)) 
listed  in  Fig.  4  causes  the  temperature -depth  profile  to  be  drawn.  The 
resultant  plot  is  shoim  in  Fig.  6.  The  surface  salinity  (which  in  this 
case  is  the  salinity  throughout  the  profile)  is  printed  to  the  left  of 
the  temperature  profile  because  the  salinity  plot  has  been  suppressed. 
Tne  generated  velocity-depth  profile  is  plotted  in  Fig.  7. 
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Example  2  uses  the  temperature  and  salinity  profile  presented  in 
Fig.  2.  The  input  deck  (listed  in  Fig.  8)  generates  the  computer  print¬ 
out  shown  in  Fig.  9  and  calls  for  four  plots:  a. temperature  profile,  a 
salinity  profile,  a  velocity  profile,  and  a  ray  trace  (Figs.  10,  11,  12 
and  13  respectively). 

SUMMARY 

The  fundamental  CONGRATS  programs,  S0990  and  S0991>  have  been 
revised  to  accommodate  temperature  and  salinity  data,  in  addition  to 
velocity  data,  as  a  function  of  depth.  Velocity  data  is  still  an 
acceptable  input,  leaving  the  programs  entirely  compatible  with  old 
input  decks.  The  ray  plotting,  eigenray  generation  and  eigenray 
processing  functions  of  the  CONGRATS  series  (as  described  in  refer¬ 
ences  (a)  and  (b))  have  not  been  affected  by  the  addition  of  the  new 
input  set.  Subroutine  INHJT  has  been  changed  in  order  to  read  and 
interpret  the  THERJ«1AL  cards  and  a  new  subroutine,  BT,  has  been  written 
to  convert  the  THERMAL  PROFILE  into  a  VELOCITY  PROFILE.  These  two 
routines  are  listed  in  the  appendix. 

JEFFREY  S.  COHEN 
Mathematician 


!ffiEI.LA  A.  GARRETT 
Math  Aid 
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TABLE  1 


FORMULA  FOR  TOE  CALCULATION  OF  SOUND  SPEED  IN  SEA  WATER 
COMPLETE  V  =  Vo  +  Va  +  Vb  +  Vc  +  Vd 
BASIC  V  =  Vo  +  Va 


In  which 


Term 


CondltionB  for 
Use  of  Term 


1493  +  3(T  -  10)  -  6  X  10'3(t  -  10)^ 

Always 

-  4  X  10"2(T  -  18)^  +  1.2(S  -  35) 

-  10“2(t  -  18)(S  -  35)  +  Z/61 

+  10"1d2  +  2  X  10"^d2(T  -  18)2  ^  io-1d0/9O 

Always 

2.6  X  10“^T(T  -  5) (T  -  25) 

T  >  25°C 

-  10"3d2(i)  .  4)(d  _  8) 

Z  >  7000  m 

1.5  X  10’3(S  -  35)2(1  -  D) 

S  <  307oo 

+  3  X  10'^I^(T  -  30)(s  -  35) 


where 

V  is  the  sound  Q)eed  in  n^s 
T  is  the  temperature  in  °C 
S  is  the  salinity  in  °/oo 
Z  is  the  depth  in  m  ,  and  D  =  Z/lOOO 
0  is  the  latitude  in  degrees 
NOTE:  Vo  can  also  be  written; 

Vo  =  1449.44  -  4.56T  -  0.046  0? 

+  1.2(S  -  35)  -  10“2(t  -  18){S  -  35)  +  Z/6l 
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TABLE  2 


Salinity  Field 
on 

Profile  Cards 

Salinity  Field 
on 

First  Card  of  Set 

Resultant  Salinity 

Zero 

Zero 

Constant  Zero 

Zero 

Positive 

Constant  Positive 

Positive 

Zero 

Profile  of  Values 

Positive 


Positive 


Profile  of  Values 
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5S.77 

55.79 


8202.0 
8561.6 
4921.2  55.91 

3926.96  56.17 

3280.8  56.33 

2624.64  56.76 

1968.48  57.06 

1640.4  57.15 

1312.32  57.24 


984.24 
656.16 
492.12 
328.08 
246.06 
164.04 
98.424 
65.616 
32.808 
0.0 


n 

[  I  I 


57.01 
57.10 
59.23 
61.16 
65.2 
65.16 
65.11 
65.21 
65.62 
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PROFILE 


FT/S 

CONTINUOUS  19.0 


39.0 


38.2 


I  il  I 


\  0 1 
1,10  il 

ill  I 


01, 

1  V 

•hi; 

or;  ^ 

»i  *  j 

'•I  . 


22 


I  i 

i,' 

'^1. 

’'l 

/il'  2 

.1 


or; 

>i'  2 


2 


‘‘t  ' 

1  r  1 
'v!  '2 

1 1';  2 


0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  1 0  0  0  '  0  0  0  0  0  0  0  0  0  0  0  0 '  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 
I  >  3  4  3  I  I  I  <  It  II  13  13  14  IS  li  nil  I)  /O  31  33  33  34  33  3i  31  31  31  30  31  33  33  34  3S  3i  33  31  3(  40  41  43  43  44  4S  «  41 II  41  SO  SI  S3  SI  S4  SS  Si  SI  Si  SS  lO  il 

11 1  Ml  I n 11 11 1 11  n  1 11 1 1  M  1 1  M  M  11111111111111111111111111111 


or;  5 

“'ll,' 


1 1 


2 


22222222222222222222222222222  22222222222222222222222  22222222 


2  2 
3 


33333  33333333,3333.33  331" 


'M  3  3  3  3  3  3  3  3  3  3  3  3  3 


Fig.  1  -  Ten^jerature  and  Salinity  Set  1 
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Fig.  4  -  Listing  of  Exan^jle  1  Run  Deck 
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Fig.  6  -  Listing  of  Example  2  Run  Deck 
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bUBKOuTiNt  BT  CONVERTS  A  TF^PhRATUKL  PROFILE  TO  A  VELOCITY 

Profile. 

Subroutine  bt(tP) 

OiMENSION  and  COM.1ON  statements. 

uImENSIOU  TARGET(1010»3) .aN6L£S(1010) ,bONARfb,2) »sUHFAC(5) » 

1  bp (210 » 2) » BOTTOM ( b) »BP(2i0.2) .TOuEKACb) »VP(210,2) » (P(210f3) » 

2  Fia(l2) 

Common  target  »  angles  »  sonar t surfac  »  sp  »  bottom . bp  »  tolera  » vp » 

1  rEVMAXiPROCES 
NCARDS  =  Tp(210»l) 

Transfer  depths  to  vp  array. 

Do  10  N=1»NCARDS 
VpCNfl)  =  TP(N,1) 

iu  Continue 

Calculate  velocity  from  temperature^  depth»  latitude»  and  salinity. 

DO  lao  Nmi.NCARDS 

Convert  depth  to  meters. 

« 

■J.  =  TP(N,1)»TP(204»1)/TP(20‘»»3) 
o  =  Z/1000 

Vo  =  1449.44  +  4.56  ♦  TP(N.2)  -  O.U4b  *  TP(N.2)  ♦  TP(N.2)  +  1.2  ♦ 

1  (lP(Nr3)  -  35.0)  -  0.01  ♦  (Tp(N»2)  -  18.0)  ♦  (TP(N,3)  -  35.0)+ 

2  Z/61.0 

VA  =  0.1  ♦  D  ♦  D  +  0,0002  ♦  D  ♦  D  ♦  (TP(N»2)  -  16,0)++2  ♦  0.1  ♦ 
ID*  TP(2i0,2)/90,0 

TcST  temperature  GREATER  THAN  OR  EQUAL  TO  25  DEG.C. 

IF  (  TP(N»2)  -  25.0  )  110,  r 

Vu  =  0,0002b  ♦  TP(N,2)  *  (TP(N»2)  -  5.0)  ♦  (TP(N,2)  -  25.0) 

GO  TO  120 
110  Vb  =  0 

TEST  DEPTH  greater  THAN  7000  M, 

120  IF  (  Z  -  7000  )  i30.  130, 

VC  =  -O.UOl  *  0  *  D  ♦  (D-4.0)  ♦  (D-0.0) 

GO  TO  140 
130  VC  =  0 

Test  salinity  less  than  30-/1000. 

140  IF  (  TP(N»3)  -  30,0  )  ,  lb0»  150 

Vo  =  0,0015  ♦  (l.O-O)  *  (TP(N,3)  "  35.0)**2  >  0, 000003  ♦  TPCN,2) 

1  ♦  TP(N,2)  ♦  (1P(N,2)  -  30.0)  ♦  (TP(N,3)  -  35,0) 

Go  TO  16u 
150  =  0 

loO  VP(N»2)  =  VO  +  VA  +  VB  +  VC  ♦  VO 
C  21 
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C  CONVERT  velocity  TO  DESIRED  UNITS, 

C 

VP(N»2)  =  VP{N»2)  ♦  TP{205,3)/VP(204»2) 

180  Continue 

Go  TO  9000 


btplot  plots  temperature  Profile  and/or  salinity  profile. 


Entry  btploT(Tp) 

IF  (  TP(208»1)  )  9000»  9000r 


a  TEST  Z  AXIS 


Compute  the  scale  factor  foR  the  depth  axis. 


call  SCAl£(TP(1»1) ,TP(208,1)»NCARDS»1»TP(207»1) ) 
1P(NCAR0S+1»1)  =  TP(NCARDS+1»1)  +  TP(NCAK0S+2» 1 )  ♦  TP(208»1) 
TP(NCARDS+2»1)  =  -TP(NCARDS+2»1) 

IF  (  TP(208»2)  }  S000»  5000»  0  TEST  T  AXIS 


PLOT  THE  depth  AXlS  IN  THE  ORIGllNAL  UNITS. 
set  the  variable  FORMAT 


FhT(1)  =  6h  DEP 
F.'1T(2)  =  6HTH  IN 
F|viT(3)  =  TP(205»1| 

call  AXlS(0.0»0.0»FMT»+lO.TP(208rl)  »90,0,TP(NCARDS.j.1»1)  » 
1  TP(NCARDS+2»1)»TP(207»1)) 

IF  (  TP(208»3)  )  ,  t  4120  ©TEST  S  AXIS 

Call  symbol(-i,o»i,o.o,i4,i2hsalinity  s  .90.0,12) 

Call  number (999,0, 999.0, 0,i4,TP(1»3) ,90.0,5) 


P1.OT  THE  temperature  AXIS  IN  THE  ORIGINAL  UNITS, 
Convert  tcmperature  to  original  units, 

COMPUTE  the  scale  FACTOR. 


4120  IF  (  TP(204»2)  -  0.75  )  »  »  4140 

Do  4130  N=1»NCAR0s 

4130  TP(N»2)  =  32,0  +  TP(Nf 2)/TP(204,2) 

4140  Call  SCAuE{TP(1,2) ,TP(208,2) ,NCARDS,1,tP(207,2) ) 


set  the  variable  format. 


FisTd)  =  6h  TEMPE 
F|vJ(2)  =  ohRATURE 
FnT(3)  =  BH  in  DE 
FmT(4)  =  6HGREES 
FhT(5)  =  Tp(2U5r2) 

Call  AXIS(0.0,TP(203,1) ,FMTr+30»TP(206,2),0.0,TP(NcARDS+l»2) , 
1  TP(NCARDS+2»2) ,TP(207»2)) 


Pi-OT  THE  TEMPERATURE  PROFILE. 


Call  line  (  TP(1,2) rTP(l»l)»NCARDS»l,l,3) 
Call  PLOT  (TP(208,2)+5,0»0.0,-3) 


HuOT  SALINITY  PROFILE. 
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C 

bOOO  IF  (  TP{208r3)  )  9000»  9000»  fi)  TtbT  S  AXIb 

Plot  depth  axis  for  salinity  profile. 

SlT  the  vARlAbLE  FORMAT 

FiviTd)  =  6h  DEP 
Fyj(2)  =  6HTH  IN 
FmT(3)  =  Tp(205»l) 

Call  axis  (0,0»0,0»FMT»+l6.Tp(20a»l)»90.0»TP(NCARDS+l»l)» 

1  TP(NCaRDS+2»1) »TP{207»1)) 

Plot  the  salinity  axis 
compute  the  scale  factor. 

call  SCALE(TP(l»3)iTP(20d»3)#NCARDS»l»TP(207.3)) 

SlT  the  VARlAbLE  FORMAT. 

FiviT(l)  =  OHSALINI 
Fiir(2)  =  6HTY  IN 
FmT(3)  =  6HPARTS 
FriT(4)  =  ohPER  TH 
F|V(T(S)  z  6H0USAN0 

Call  AXIS(0.0»TP(20S,1) »FMTr+30»TPl208,3) »0,0»TP(NcARDS+1»3) » 
1  rp(NCARDS+2»3) #TP(207»3)) 

PuOT  THE  salinity  PROFILE. 

Call  LINE  (TP(l»3i»TP(ia),NCAR0Srl»l»3) 

Call  PLOT(TPC2o8»3)+5.0»0.0»-3) 

9000  Return 
End 
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SuBKOUTINfc.  INPUT 


Dimension  tak6ET(10io»3) »aN6leS(1oio) ,soNAR(b,2) »surfac(5)» 

1  Sp(2l0»2) »B0T10M(5) »BP (210 »2) » TOLEKA ( b ) » VP (21 J , 2) » TP (210 »3) 

Dimension  Gf<oup(  h  ) /obHcoMMENTAROEi bonak  sukfacbottomvelocitherma 
IPKoCESEND  FREQUEmAXIMU/ 

dimension  c00H0(8)/4aHRANGE  DEPTH  ANGLE  LOSS  TOLERAAXES  PROFILP 
ihase  / 

dimension  units (2, 32)/ 


1 

i.0936111lE-5, 

6hCM  » 

1.0936HllE-b. 

6HCM/S 

9 

2 

2, 7777777 7E-b, 

6h1N  r 

2,77777777E-*5» 

dHIN/S 

9 

3 

u •33333333E~3» 

6hFT  » 

0 , 33333333E**3  • 

bHFT/S 

$ 

4 

1.00000000E-3» 

6hYD  » 

1.00000i)OOt-3» 

bHYD/S 

9 

5 

1.09361lll£-3, 

6hM  r 

1.093bllllE-3. 

6HM/S 

9 

6 

*i.00000000E-3» 

6riF  » 

2,OUOUOoOOE-3, 

oriF/S 

9 

7 

U.33333333E-0, 

6hKFT  t 

0,33333333E-0» 

6HKFT/S 

9 

a 

1.00000000E-0» 

6HKY0  » 

l.OOOUOOOOE-0. 

6HKYD/S 

9 

9 

1.09361111E-0, 

6hKM  f 

1.09361111E-0. 

bHKM/S 

9 

T 

i.76000000E-0» 

6hMI  0 

l»7bOOOOOOE-0» 

6HMI/S 

9 

1 

^.026B0000E-0| 

6hN  Ml  0 

2.026bOOOOE-O. 

6HKN0TS 

9 

2 

i.74532925E-2, 

6hOEG  0 

1.74b3292bE-2» 

6HDEG/S 

9 

3 

i,ooooooooE-o» 

UhRAD  0 

l.OOOOOOOOE-0. 

6HRA(J/S 

9 

4 

b.28318b3lE-0» 

6iiCPS  0 

6,2b3ia53lE+3» 

bHKCpS 

9 

5 

l.UOOOOOOoE-3, 

6riMS  0 

l.OOOOOOOOE-0. 

6HSEC 

9 

6 

1.00000000E-0» 

bHC  0 

5,55bb5fi5bE"l» 

6HFAHR 

/ 

DIMENSION  TEST(l2)»FMT(l2)»DATA(204b) ,REVMAX(4) 


UlMENSION  CMMNT(12)/72H  THERE  SHOULD  QE  AT  LLAST  dNE  COMMENT  SET 
1.  / 
Data  PRINT/6hPRlNT  / »LlNMAX/60/f ICMMNT/0/ 

Common  target  »  angles i sonar  »  surfac  r  sp  »  bottom  »  dp . t  olera . vp » 

1  hEVMAX»PtS0CES 


set  Initial  conditions. 


A,^Ol£S(1010)  =  0.0 
TARGETdOlOfl)  =  0.0 
rARGEr(ioio»2)  =  0.0 

IhEAD  s  0 

LINES  =  lInMAX 

NhRINT  =  FlD(33»3,PROCES) 

NiAPEl  =  FLD(30r3,PR0CES) 

NfAP£2  =  FlD(27»3,PR0C£S) 


Check  the  group  code. 


100  Read  loi,  (TEbT(J) » J=i»9) 

101  FoRMAT(  3(A6»A4)»  3F10.5  ) 

DO  110  I=lrll 

IF(  TEST(1).GT.GRoUP(I)  .or.  TEST(1),LT.6R0UP(I) 
loROUp  =  I 


)  oO  TO  110 
24 
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Go  TO  (lQU0»200»2o0#200»2O0»20U»200,a000f9000»300»g200) »  IGROUP 

110  Continue 

Th£  group  cODc  is  INCORRECT. 

120  KiviTd)  =  (+6HGR0UP  ) 
iJO  PkInT  132»  FMld) 

1^2  format (  lOX  IAHTHe  following  ,  A6,  ItiHcOOt  IS  INCORRECT.  ) 

140  Print  142»  (  TESTiJ) .J=1»9  ) 

142  FoRMAT(  24X  6A6»  ^Fl2.5  ) 

ISO  print  152 

-152  Format (  lOx  26hthe  program  cannot  continue.  ) 

Stop  6 


Check  the  coord  code. 


200  Go  210  I=lf8 

IK  T£ST(:5).G7,COoRO(1)  .OR.  TEST  (S)  .LT. COORD  ( I )  )  GO  TO  210 
ICOORO  =  1 

Gu  TO  (300. 300, 500.500. 300. 3oo»600»b00) r  iCOORO 
210  CONTINUE 

The  COORD  CODE  IS  INCORRECT. 

F|.Jd)  =  (+6HC00RD  ) 

GO  TO  130 


Check  the  units  code. 


300  Do  310  1=1.30 

IF(  TEST(5),GT,UNiTS(2.I)  .OR,  TEST(5) ,LT.UniTS(2. I )  )  GO  TO  310 
IjNlTS  =  I 

1F(  IGROUP.EQ.IO  )  60  TO  320 

1F(  ICOORO, GE. 6  )  GO  TO  400 
1F(  TEST(9}  )  305.315.305 

305  hoards  =  (TESTC8)-TEST{7))/TEST(9)  +  1.5 
1F(  NCARD5*(1001-;mCARUS)  )  2100.2100.320 

310  Continue 

TnE  units  CODE  IS  INCORRECT. 

FnTd)  =  t+oHUNlTs  ) 

GO  TO  130 

Print  the  data, 

315  TLST(a)  =  TESK?) 
imCARDS  =  1 

320  iF(  LINES. LT.LINMaX-S  )  GO  TO  330 
PrvINT  322 
322  FuRMAT(  IHI  ) 

Lines  =  i 
GO  TO  340 

330  1F(  IHEAO.EO.I  )  GO  TO  350 

Print  332 
332  format (  //  ) 

LINES  =  lines  4  3 

340  Print  342 

342  Format  (  2i>x  73HGROUP  COORDINATE  UNITS  INITIAL 

1  final  Increment.  /  )  25 
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LINES  =  lInES  +  2 
IhEAO  =  1 

slt  The  vaf^able  format. 

OOU  FmT(1)  =  (+6H(  24x  ) 

P|^T(2)  =  (+6HbAfa»  3) 

IF(  UiMITS(1»IUNITs)  .GT.2.0E-3  )  GO  TO  360 
FmT(3)  =  (+6HF12.2)) 

60  TO  370 

36U  FmT(3)  =  (+6HFl2,b)) 

370  PkINT  FMT»  (TtST(J} »J=1»9) 

LINES  =  lines  +  1 

TLST(7)  =  TEST(7)  ♦  uNlTS(l»IuraTS) 

TtSTtd)  =  TEST(8)  ♦  UNITS(1»IUNITS) 

(uST(9)  =  tES7(9)  *  UNITS(1»IUNITS) 

6u  TO  (3d0»2000»3o00. 4000, 5000»6000»36o»360»:^b0»9100»  100) »  I6R0UP 

The  codes  are  inconsistent. 

3bU  PkINT  382 

3d2  Format {  lUx  33hthe  above  codes  are  inconsistent.  ) 

60  TO  150 


Check  the  axes  card. 


400  IF(  NTAPEI.EQ.O  )  NTAPElsl 
iF(  IUNITS.GT.4  )  GO  TO  420 
ipAPER  =  600  -  50o  *  { (IUNITS-1)/2) 

Uu  4l0  «Jzl,3 

1F(  TEST(d+6)  «  (i00.0*UNITS{l,3)-T£ST(J+b)*UNITS(i,IUNITS) )  ) 
1  420^410, 4l0 

410  CONTINUE 

6o  TO  (420, 420, 420, 4100, 5100,6100, 7100, 420r420) ,  IGROUP 

TmE  axes  Card  is  incorrect. 

420  PkINT  422 

422  Format (  lox  37hthe  following  axes  card  is  incorrect.  ) 

60  TO  140 


RcAD  THE  LOSS  TABlE. 


500  IhEAU  =  0 

LINES  =  lInES  +  2ti 

1F(  LINES, lT.LINMaX-S  )  60  TO  510 

PkINT  322 

LINES  =  25 

6u  TO  520 

Dili  PkINT  332 

520  IF(  IGROuP.EQ.4  )  GO  TO  540 
IF(  1GR0UP.NE,5  )  GO  TO  120 
1F(  IC00RD.EQ,4  )  GO  TO  530 
TcSTCl2)  =  BPHASE(FREQ) 

60  TO  100 

530  TLSr(12)  =  BLOSS(FREQ) 

6o  TO  100 

540  IF(  IC00RD.EQ,4  )  GO  TO  550 
TlST(12)  =  SPHaSE(FREQ) 
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Go  TO. 100 

550  TcST(12)  =  SL0SS(|.REQ) 

Go  TO  100 


Check  the  numeer  of  cards  in  the  PKOFiLt. 


600  NCAROS  =  TEST(7) 

IK  (NCAROs.6t.2).AND.(200.GE,NCARUS)  )  GO  fO  610 
PkInT  602»  (TEST(j) ij=l,4) ,  NCARDS 
602  FuRMAT(  iUX  27HTHE  NUMBER  oF  CARDS  IN  THE  t  2(A6rA4)» 

1  110»  44H»  exceeds  200  CARDS  OH  IS  LLbS  THAN  2  CAKUS.  ) 
Go  TO  150 

PhINT  THE  heading. 

6i0  InEAO  =  0 

lines  =  ncaRDS/2  +  10  +  lines 
IK  LINES. uT.LINHaX-1  )  GO  TO  620 
PkINT  322 

UNES  =  nCaRDS/2  +  6 
Gu  TO  630 
02O  print  332 

630  print  632#  (TEST! J) » J=l,4) 

632  PoRMATl  5iX  2(A6»a4),  /»  5lX  17H - - - -  /  ) 


Check  the  units  code. 

700  HtAD  101,  TEST(l),TEST(l0),TEST(2)»TEST(ii)»TEST(3),TEST(12) 

J  =  I 

710  Do  730  1=1,32 

1P(  T£ST(U) .GT,UNlTS(2,I)  .OR,  TEST  10) .LT.UNITS(2» I)  )  GO  TO  730 
IF(  J.GT.l  )  GO  10  720 
luNiTS  =  I 
0  =  2 
Go  TO  710 

720  IF  U  ,GT.  2)  GO  TO  725 
OuNlTS  =  I 
0  =  3 

IF  (XGROUP  .EQ,  7)  GO  TO  7IO 
Go  TO  740 
725  RuNITS  =  1 
Go  TO  740 
730  Continue 

c  TmE  Units  code  is  incorrect. 

P(vJ(l)  =  (+6HUNITS  ) 

Pr^lNT  132»  FMT(1) 

Print  14*;»  TEST(1),TEST{3),TEST(2),TEST(4) 

'  GO  TO  150 

c 

c  StiT  THE  variable  FORMAT. 

740  FiiTd)  =  (+6H(  2(  ) 

FmT(2)  =  (+6H  115,  ) 

IF(  UnITs(1»IUNITs).GT.2.0E-3  )  GO  TO  750 
FmT(3)  =  (+6HF14.2,) 

Go  TO  760 

750  P,.,T(3)  =  (♦6HF14.5#) 
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760  if-  (IGROUP  •E'i.  7)  GO  TO  763  . 

LuNiTS  =  JUNITS 
Go  TO  765 

7b3  LUNITS  =  KuNITS 

7ob  (  UNIT6(1»LUN1tS)  .GT.  2.0£i:-3  )  Gu  .VO  77(1 
f-HT(4)  =  (+6HF16.2») 

Go  TO  780 

770  FmT(4)  =  (♦6HF16,S») 

780  F|mT(5)  =  (+6H15X)  )) 

GO  TO  (380,380»3a0»4200»52Q0,6200»72U0.380»3rt0) »  IGROUP 


HlAO  and  WRlTt  A  COMMENT  CARD. 

1000  KtAO  100l»  (TLST(J)»0=1»12) 

1001  FlRMAT(  i2A6  ) 

1002  FoRMAT(  iOx  12A6  ) 

1F(  lines. LT.LINMaX  )  GO  T®  1005 
HkINT  322 
LINES  =  1 

100b  HiUNT  lOO^t  (1EST(J)«J=1»12) 

LiNES  =  lines  +  1 
loMMNT  =  IcMMNT  +  I. 

1F(  ICMMnT.GT.1  )  GO  TO  lOO 
Uo  lOlO  o=lfl2 
1010  C|.iMNT(J)  =  TEST(J) 

Go  TO  100 


Convert  and  store  the  target  data, 

2000  IF(  IC00RD,GT,2  )  GO  TO  380 
1F(  NTAPE2,EQ,0  )  NTAPE2=1 
TmRG£T(1004»IC00R0)  =  UNlTSdflUNITS) 
(MRGEr(l005»IC0ORu)  =  UN1TS<2» lONlTS) 

N  =  target ( 1010 rIcOORD) 
iMbAROS  =  N  +  UCAHOS 
lr(  NCARoS.GT.lOOo  )  GO  TO  2100 
TARLET(l0i0»lLOORu)  =  NCARDS 
=  N  +  1 

UU  2010  i=n»ncaros 

2010  TmRGEKI.IcOOKO)  =  TEST{7)  ♦  ( X-N) ♦TEST (9) 
GO  TO  100 

The  number  of  entries  has  been  exceeded. 

2100  PkINT  21uii,  GF^OUP  ( IGROUP)  r  COOKD  ( ICOORO) 

2102  Format (  lox  i^hthe  number  of  »  a6»  ixa5» 

'1  <.3HS  exceeds  lOOo  EFmTRIES,  ) 

Go  TO  150 


Convert  aI'ID  store  the  sonar  data. 


5000  XF(  ICOORD-3  )  3l00.3200»380 

5100  SoNAR(lrlCOORD)  =  TEST(7) 
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C 

C 

3200 


3210 


C 

c 

3300 

3302 


C 

c 

c 

c 

c 

^000 


c 

c 

4100 


c 

c 

c 

c 

4200 


4202 

4212 


SONAR(2»ICOORO) 
SyNAR(3»iC00RU) 
S0NAR(4f ICOORD) 
SoNAR(5»ICOORD) 
GO  TO  100 


=  T£ST(8) 

=  TEST(9) 

=  UNITS(1»IUNITS) 
=  UN1TS(2»IUN1TS) 


Si  ORE  THE  SONAR  AnGUE  DATA, 

1F(  AcJS(Tfc.ST(7) )  .gT. 1.5707  .OH.  TEST(6)  .GT. 1.5707  ) 
AuGLES(1004)  =  UNiTS(1»IUNiTS) 

Angles (1005)  =  unitS(2»iunits) 

N  =  ANGLE5(1010) 

NCAROS  =  N  +  NCAROS 

XF(  NCARQS.GT.IOOo  )  GO  TO  2100 

A,s.GLES(1o10)  =  NCaROS 

N  =  N  4-  1 

Do  3210  I=N»NCARqS 

AnGLES(I)  =  TEST (7)  +  ( I-N)*TEST (9) 

Go  TO  100 


GO  TO  3300 


The  SONAR  angles  exceed  their  BOUNDS. 

HhINT  3302 

Format {I0X5OHALL  sonar  angles  must  lie  between 

Go  TO  150 


-90  AND  4-90  DEG,} 


Convert  and  store  the  surface  data. 

IF(  ICOORD, NE, 2  )  GO  TO  380 
SuRFAC(I)  =  TEST (7) 

SURFAC(2)  =  TEST(a) 

SuRFAC(3)  =  T£ST(9) 

SuRFAC(4)  =  UNiTSdrIUNiTS) 

SuRFAC(5)  =  UNITS(2dUNlTS) 
bo  TO  100 

STORE  the  surface  AXES  DATA. 

SH(207»1)  =  2.777777777E-4/UNlTS(l»iUNlTS) 
Sh(207»2)  =  SP(207»1) 

Sp(20a,l)  =  TEST(7) 

Sp(203»2)  =  TESTCa) 

Go  TO  5100 


rlad  and  Print  the  surface  profile, 

Sp{204,l)  =  UNiTSd.lUNiTS) 

Sh(205»1)  =  UNITS(2»IUNITS) 

Sp(204f2)  =  UNlTSd»jUNlTS) 

SF(205,2)  =  UNITS(2»JUNITS) 

Sp(2l0»l)  =  NCARDS 
Sh(210»2)  =  TEST(5) 

SuKFAC(4)  =  0.0 

FnINT  4202»  SP(205»1) »SP(205»2) »SP(205,1) »SPi205»2> 

FoRMATC  2(12X  IGHCARD  RANGE-»  A6»  9H  DEPTH-.  Ab.  11X)»  /  ) 

KlAD  4212»  (SP(N»1) rSP(N»2) »N=1»NCARDS) 

FuRHAT(  2F10,5  ) 

J  =  (  NCaRDS  +  1  )  /  2 
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UO  4220  U=lfJ 
(><  =  L  +  J 

IF(  N.LE.NCARDS  )  GO  TO  4220 
PkINT  FMT»  L»SPCL,1) ,SP(L»2) 

Go  TO  4230 

4220  PkINT  FMT»  L»Sp(L. 1) ,SP(L»2) »  N»SPtN» 1 ) »bP (N.2) 

C 

C  CONVERT  and  CHECK  THE  PROFILE, 

4230  Do  4260  N=l»NCARoS 

SP(N»1)  =  SP(N,1)  ♦  SP(204fl) 
bP(N»2)  =  SP(N»2)  ♦  SP(204,2) 

C 

c  lb  the  range  increasing. 

IF(  N.EQ.I  )  GO  TO  4260 
iF(  SP(N,1).G1.SP(N-1»1)  )  GO  TO  4260 
NCAROS  =  N 
Go  TO  427O 

4260  CONTINUE 
GO  TO  100 
C 

C  A  CARD  Ib  INCORRECT  OR  OUT  OF  ORDER. 

4270  print  4272,  NCARDS 

4272  Format (  lOX  IIHCArD  number,  I4»  13hIS  INCORRECT.  ) 

Go  TO  150 
C 
C 
C 

C  CONVERT  and  STORE  THE  BOTTOM  DATA, 

Q  - - - - - - 

6000  IF(  IC00r0.NE,2  )  GO  TO  380 
boTTOM(l)  =  TEST(7) 
boTTOM(2)  =  TEST (a) 
buTTOMO)  =  TEST  T9) 

BoTTOM(4)  =  UIMITS(1,IUNITS) 
boTTOM(5)  =  UNITS(2»IUNITS) 

Go  TO  100 
C 

C  blORE  THE  BOTTOM  AXES  DATA. 

6100  bP(207»l)  =  2.777777777E-4/UNITS(l,iUNlTS) 

bp(207,2)  =  BP(207il) 
bp(20a,l)  =  TEST(7) 

Dp(208,2)  z  TESTta) 

Go  TO  100 
C 
C 

c  Read  and  print  the  bottom  profile, 

Q  - - - 

o20U  bp(204,l)  z  UNITS(1»IUNITS) 
bp(205,l)  z  UfUTSc2»lUNlTS) 

'  bp(204,2}  z  UniTS(1»JUNITS) 

bp(205,2)  z  UNITS{2»JUNITS) 
bR(2l0,l)  z  NCARDS 
op(2l0»2)  z  TESTIS) 
buTrOM{4)  z  0.0 

Print  42u2,  BP(20b»l) »BP(2o5,2) »DPC205,1) »BPf206»2) 

KeAD  42l2r  (BP(N»l) »bP(N»2) »Nzl»NCARDbJ 
J  z  (  NCaRdS  ♦  I  )  /  2 
DO  5220  L=1»J 
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N  =  L  +  j 

IF(  N.LE.WCAROS  )  GO  TO  5220 
HkINT  FMT»  L»BP(L,1) ,bP(L»2) 

Go  TO  5230 

b22U  P.atMT  FMTr  U»uP(L,  1)  »BP(L,2) »  N»t3PlN,l>  ,tiP(N,2) 

C 

c  Convert  and  check  the  profile. 

5230  Do  b260  N=1»I'jCARdS 

Dp(N»l)  =  BP(N»1)  ♦  aP(204.1) 
bp(N»2)  =  BP(N»2)  *  8P(204»2) 

C 

C  15  the  range  IfjCREASING, 

1F(  N.EQ.I  )  GO  TO  5260 
IF(  aP(N,l) .GT,BP(N-1»1)  )  GO  TO  5260 
NCARDS  =  N 
Go  TO  4270 
5260  CONTINUE 
Go  TO  100 
C 
C 

c 

c  Convert  and  store  the  velocity  data, 

Q  - - — - - - — 

oOOO  IF(  ICOOrD.NE.S  )  GO  TO  380 
TOLERA(l)  =  TEST(7) 

ToLERA(2)  =  TEST(a) 

ToLERA(3}  =  TtST(9) 

ToLERA(4)  =  UNITS(1»IUNITS) 

TOLERAO)  =  UMTS(2»IUNITS) 

Go  TO  100 
C 

C  STORE  THE  VELOCITY  AXES  DATA, 

biuU  VP(207,1J  =  2,777777777E-4/UNITS(1»IUNITS) 

VF(207»2)  =  VP(207»1) 

VP(206,1)  =  TEST(7) 

Vp(20a»2)  =  TLST(a) 

Go  TO  100 
C 
C 

C  KlAD  and  print  the  velocity  PROFILE, 

Q  - - - - - - - - 

6200  Vp(204»l)  =  UUITSdflUNlTS) 

Vp(205»l)  =  UNITS(2» IUNITS) 

Vp(204,2)  =  UNITS(1»JUNITS) 
vp(205,2)  =  UNITS(2»JUNITS) 

Vp(2l0d)  =  NCARDS 
Vp(2l0»2)  =  TLST(5) 

Print  6202,  VP(205»1) fVP(205,2) »VP(205,1) »vp(205»2) 

6202  FoHMAT(2(l2X  16HCAR0  DEPTH-r  A6»  lOH  VELOCITY-.!  A6*  lOX)  i  /) 

'  read  4212»  (VP(N»1) »vP(N»2) »N=1»NCARDS) 

621U  J  =  (  NCaKuS  +  1  )/2 
DO  6220  LzlfJ 
N  =  L  +  J 

IF(  N.LE. NCARDS  )  GO  TO  6220 
PRINT  FMT»  L»VP(L,1) ,VP(L»2) 

Go  TO  6230 

6220  PrxINT  FMT»  L»VP(L,1)  ,  VP{L,2)  #  N!VP(N»l) ,  VP(N,2) 

C 
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c  Convert  and  check  the  profile, 

6230  Do  6260  N=lfNCARQS 

Vp(N»l)  =  VP(N,1)  4c  VP(204.1) 

Vp(N»2)  =  vP(N»2)  ♦  wP(204»2) 

IF(  VP(N.2)  .GT.l.bS  .AND.  1.77,6T, VP(N.-2)  )  00  TO  6240 

TtiE  velocity  exceeds  ITS  BOUNDS. 

PniNT  6232»  N 

6232  format {  lOX  28hVEl0CITY  PROFILE  CARD  NUMBER#  14, 

1  tLOti  exceeds  its  rjOUNOS,  ) 

00  TO  ISO 

IS  the  depth  increasing, 

6240  1F(  N,EQ.1  )  GO  TO  6260 

IF(  VP(N,1) .GT,VP(N-1,1)  )  GO  TO  6260 
NcARDS  =  N 
Do  TO  4270 

6260  Continue 
Go  TO  100 


Convert  aNo  store  temperature  data. 


STORE  THE  THERMAL  AXES  DATA. 

Tp(207,l)  =  2.777777777E-4/UNITS(1,1UNITS) 
Tp(207,2)  =  TP{207,1) 

Tp(207r3)  =  TP(207,1) 
lp(208,I)  =  TEST  (7) 

TP(208,2}  =  TEST  (8) 

Tp(20d,3)  =  TEST  (9) 

GO  TO  100 

read  and  Print  the  thermal  profile. 


720u  vp(204,l)  = 

VP(206,1)  = 
Vp(204r2)  = 
Vp{205,2)  = 
Vp(2l0,l)  = 
Vp(2l0»2)  = 
TP(204,1)  = 
lP(205,l)  = 
Tp(204,2)  z 
rpi205,2)  z 
,  Ip(204,3)  z 

Tp{20S,3)  z 
TP(210rl)  z 
Tp(2i0,2)  z 
TP(210,3)  z 
PkINT  7202, 


UNITS(1,1UNITS) 

UnITS(2,IUNITS) 

UNITS(1,KUNITS) 

UI^ITS(2,KUNITS) 

NCARDS 

TEST(5) 

UfaTS(l,lUNlTS) 

UNITS(2,IUNITS) 

UNITS(1,JUNITS) 

U.NITS(2,JUNITS) 

UnITS(1,9) 

UNITS(1,10) 

NCARDS 

TEST(6) 


TEST  (9) 

TP(20b,l),TP(205,2),TP(20b.l),TPj205,2) 

7202  Format  ( 2  (6X,16HCAKD  DEPTH-,A6,12HTEMPERATURE-.Ab,14HSAHNITK 

1-/1000, 6X),/) 

RcAD  720b»  (TP(N,i) ,TP(N,2)iTPJN, 3) ,Nzi, NCARDS) 

720b  FoRMATOFlO.b) 
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Tt;sT  salinity  input  to  use, 

SaL  =  0,0 

Uo  7207  1  =  1»NCARDS 

IF  (  TP(1.3)  )  7208#  f  f 

SAL  =  sal  +  TP(I»i) 

If  (SAL  .GT.  0,0)  GO  TO  72l5 

7207  Continue 
Go  to  7210 

7208  PhINT  7257,  I 
print  7209 

7209  Format  (10X»36MNEgATIVE  salinity  is  not  valid  data.) 

Go  to  150 

7210  uo  7212  I  =  1»NCArDS 
Tp(I»3)  =  TP(210»3) 

7212  Continue 

Print  the  temperature  profile, 

721b  J  =  (  NCaRDS  +  1  )/2 

Uo  7220  L=lrJ 

N  =  L  +  J 

IF  (N  ,LEt  NCARDS)  GO  TO  7220 
Print  722^.  l»tp(l»1) »tp(l»2),tP(l»3) 

GO  TO  7226 

7220  Print  7225,  L»TP(L»l)»TP(L»2),TP(L»3)»N»TP(tJ,i)»TP(N»2)»TP{N»3) 
722b  FoRMAT(2(bXrI3,6X,F8.2»aX,F7,2»llXrF7,2»10X) ) 

7220  Print  7227,  TP(21o»2) 

7227  FoRMAT(//»45X»12HlATITUDE  =  »Fl0,b»bH  DEGREES,//) 

CONVERT  and  CHECK  TEMPERATURE  PROFILE 

7230  Do  7260  N=lfNCARDs 

Ip  (wlUNITS  ,EQ,  3i)  GO  TO  7235 
TP(N#2)  =  (TP(N»2)  -  32.)  ♦  TP(204r2) 

723b  IF  (TP(N,2)  ,GT.-3,0  ,AND.  35,0  .GT,  TP(N»2))  GO  TO  7250 

IhE  TEMPERATURE  ExCEEDS  NEpTUNlAN  bOUNoS, 

Print  7237,  n 

7237  Format (lOX,31hTEMpERATURE  PROFILE  CARD  NUMBER » 14 » 19HEXCEE0S  ITS  B 
IOUimDS.) 

check  if  salinity  IS  within  NEPTUNIAN  BOUNDS. 

72j0  If  (SAL  .GT.  0,0)  Go  TO  7254 
K  =  N 

IF  (K  ,GT.  1  )  Go  TO  7260 

iP  (TP(2iO,3)  ,GE.  0.0  .ANDf  43,0  ,GT.  TP(210»3))  GO  TO  7260 
Go  TO  7256 

7254  IF  (TP(N,3)  .GE,  0.0  .AND,  ‘♦3,0  .GT.  Tp(210,3))  Go  TO  7260 
ThE  salinity  EXCEEDS  NEPTUNIAN  BOUNDS. 

7255  Print  7257 

Print  7237,  n 

Go  TO  7260 
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72bfa 

72by 

72oJ 


7320 

7330 

733b 


aooo 


<J0iO 

b02U 

c 

c 


c 

c 

tiluo 


c 

c 

0200 


0210 

0220 

0230 


0240 
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PKINT  72b7 

FoRMAT(10X,40hSALlNnY  IS  NOT  WlfHiN  NEPTUNIAN  BOUnOS.) 

Continue 
Call  bt(tP) 

pkint  thc  heading 

IhEAO  =  0 

lines  =  nCaRds/2  +  10  +  lines 
Ip  (LINES  .LT.  LInMAX  -  1)  60  TO  7320 

print  322 

LINES  =  nCaRDS/2  +  6 
GO  TO  7330 
Print  w 
Print  7335 

Format  (  3ix»17Hvelocity  profile  »/»51x»i7h - — - - »/) 

print  the  velocity  profile. 

Print  0202,  VP(205»1) »VP(2o5,2) »VP(20b.l) »VP(205»2) 

Go  to  6210 


Process  congrats. 


IF(  PRINT. lT.TEST(3)  ,0R.  PRINT. GT.TEST(3)  )  GO  To  6010 
NpRiNT  =  1 
GO  TO  8020 
NpRINT  =  0 

KLD(27f9,PR0CES)  =  NpRiNT  ♦  8  *  (NTAPEl  +  8*nTAPE2) 

IS  there  At  least  one  comment  card. 

IK(  ICMMNT.GT.O  )  GO  TO  8l00 
Print  loui.  (CMMNr(J) »J=i.i2) 
lines  =  lInes  +  1 

liJTlALl2E  THE  PLOTTER  TAPE. 

IF(  NTAPEl. NE.l  )  GO  TO  8200 
C«LL  PL01S(  data, 2048.1  ) 
call  PLOT(0,(j»0,0»  IPAPER) 

Call  ploT(5,o»o.o»-3) 

CmLL  symbol (0,0. U.O, 0.14, cMMNT (1) »9U,0» 72) 
call  PLOT{5,0.0,0.-'3) 

UlAPEl  =  2 

Sort  the  Target  coordinates. 

Du  8260  j=1.2 

N  =  target (1010. J) 

IK(  N.LT,2  )  60  to  8260 

L  =  1 
1  =  L  +  1 

IF(  TAR6ET(L.J)-TaRGET(I.J>  )  6260.8230.8240 

TaRGET(I,J)  =  TAR6ET(N.J) 

(4  =  N  -  1 
Gu  TO  8250 

T,,RGET(N+1.J)  =  TaRGET(I»J) 

TaRGET(I.O)  =  TARoET(L.J)  34 
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TaR6ET(L,J)  =  TARgET(N4-1»J) 

8260  1  =  I  1 

IF(  I.LE.lM  )  GO  TO  8220 
L  =  L  ^  1 

IF(  L.LT.N  )  GO  TO  8210 
lARGET(lOiO»J)  =  N 
a2ou  Continue 

call  UTPlOtCTP) 

c 

c  slt  the  maximum  reversal  increment.  \ 

1F(  REVMAA(2) ,GT.kEVMAX(1)  )  GO  TO  8270 
Kt.VMAX(3)  =0.0 
Go  TO  6300 

8270  N  =  TARGET (10l0»X) 

Kc-VMAX{3)  =  (rARGF.T(Nrl)-TARGE:T(i»l))/fREVMAX(2)-ReVMAX(l)) 
C 

G  SuRT  THE  Sonar  angles. 

0300  U  =  AN6LES(10i0) 

ih(  N.LT.I  )  go  to  8400 
Do  8305  L=1»N 

iF(  A8S(aNGLES(L) ) .GT.I.OE-4  )  GO  TO  8305 

Angles (L)  =  o.o 

8306  Continue 

IF(  N,LT,2  )  GO  TO  8400 
L  =  1 

8310  I  =  L  +  1 

8320  IF{  angles (L)-ANGlES( I)  )  8350 rG330 »8340 

83 jO  AnGL£S(I)  =  ANGLES (N) 

N  =  N  -  1 
Go  TO  8350 

8340  AijGLES(N+l)  =  ANGlES(I) 

augles(I)  =  angles (L) 

Angles  (L)  =  angles  (N+d 

8350  1=1+1 

iF(  I.LE.N  )  GO  TO  8320 
L  =  L  +  1 

1F(  L.LT.N  )  GO  TO  8310 

Angles (1010)  =  n 
8400  return 
c 
c 
c 

c  Terminate  congrats. 

c  - - 

9000  Print  9002 

90 0*;  Format  (  1H1»  9X  28HC0N0RATS  HAS  BeEN  COMPLETED.  »  /»  iHl  ) 
9030  iF(  NTAPE2,EQ.0  )  GO  TO  9090 

End  file  2 

9090  S|0P  5 

c 

c 

c 

c  READ  miscellaneous  DATA. 

^  - - - 

C 

c  read  THE  Frequency. 

9100  TaRGeT(1008»3)  =  tEST(7) 

IF(  (IUNITs-2/)+(iUNITS-28)  )  380»9110»380  35 
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'■JixU  rAKGET<lu08»l)  =  ATTtN(TAKGfc.T(100a»5)) 

PkINT  91l2f  TARGETdOOO.l) 

9U2  FoRMAT(  12HA  iTErjUATlGN  t  1  bHUb/KYUf  oX  f-l2.b  ) 

LINES  =  lUJES  ^  1 
Go  TO  100 
c 

C  HEAD  THE  MaXTK.uM  NUMBER  OF  REVERSALS, 

9200  RdVHAX(l)  =  TEST (7 ) 

KEVMAX(2)  =  MAX(TEST(7).TEST(6)) 

RlVMAX(4)  =  TEST(9) 

Tt:ST(8)  =  REVMAX(2) 
lufJlTS  =  lb 

Xh(  REVMaX(2) ,LT. 201.0  )  GO  TO  320 
Print  9202 

9202  FoRMAT(  l0x45HTHE  MAXIMUM  NUMBER  OF  KEoERSALS  EXCEEDS  200.0  ) 
Go  TO  140  ! 

fnd 


